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We have studied big bang nucleosynthesis in the presence of regions of antimatter. Depending on 
the distance scale of the antimatter region, and thus the epoch of their annihilation, the amount of 
antimatter in the early universe is constrained by the observed abundances. Small regions, which 
annihilate after weak freezeout but before nucleosynthesis, lead to a reduction in the 4 He yield, 
because of neutron annihilation. Large regions, which annihilate after nucleosynthesis, lead to an 
increased 3 He yield. Deuterium production is also affected but not as much. The three most impor- 
tant production mechanisms of 3 He are (1) photodisintegration of 4 He by the annihilation radiation, 
(2) p 4 He annihilation, and (3) n 4 He annihilation by "secondary" antineutrons produced in 4 He an- 
nihilation. Although p 4 He annihilation produces more 3 He than the secondary n 4 He annihilation, 
the products of the latter survive later annihilation much better, since they are distributed further 
away from the annihilation zone. Our results are in qualitative agreement with similar work by 
Rehm and Jedamzik, but we get a larger 3 He yield. 

PACS numbers: 26.35.+C, 98.80. Ft, 98.80.Cq, 25.43. +t 



I. INTRODUCTION 
A. Antimatter in the Universe 

The local universe is baryon asymmetric. It contains 
matter, not antimatter. In standard homogeneous big 
bang cosmology the universe was filled with a uniform 
mixture of antimatter and matter very early on, with a 
slight excess of matter over antimatter. This excess of 
matter was left over when matter and antimatter anni- 
hilated during the first millisecond. 

This baryon asymmetry is characterized by the baryon- 
to-photon ratio, 



where rib is the number density of baryons and the 
number density of antibaryons. 

There are many proposed mechanisms for baryogenesis 
to explain the origin of this asymmetry. The simplest ver- 
sions of baryogenesis produce a homogeneous asymme- 
try, but there are many possibilities for inhomogeneous 
baryogenesis, which could produce a baryon excess in 
some regions and an antibaryon excess in other regions. 
This leads to a structure of matter and antimatter regions 
after local annihilation during the first millisecond. 

In scenarios connected with inflation, there is no a pri- 
ori constraint on the distance scale of these matter and 
antimatter regions. If the distance scale is small, the an- 
timatter regions would have annihilated in the early uni- 
verse, and the presence of matter today requires asym- 
metric baryogenesis, producing more baryons than an- 
tibaryons. If the distance scale is large, antimatter re- 
gions would have survived till present. 



In the latter case, an overall baryon symmetry re- 
mains a possibility. The universe could contain equal 
amounts of matter and antimatter, spatially separated 
into matter and antimatter domains. In this case, the 
absence of observed annihilation radiation from the do- 
main boundaries indicates that the typical size of these 
domains would have to be very large. 

Considering only the conditions in the present uni- 
verse, the lower limit to the domain size corresponds to 
the scale of cluster of galaxies, of the order of 20 Mpc [0. 
Because of the low density of intergalactic space between 
clusters the annihilation radiation between a cluster and 
an "anticluster" could have escaped detection. 

However, the isotropy of the cosmic microwave back- 
ground (CMB) rules out large voids between matter and 
antimatter regions during an earlier time. Thus anni- 
hilation would have been more intense before structure 
formation. The relic 7 rays would contribute to the cos- 
mic diffuse gamma (CDG) spectrum. The observed CDG 
spectrum gives a much larger lower limit to the domain 
size, of the order of 10 3 Mpc, comparable to the size of 
the visible universe |^| . A boundary of an even larger do- 
main intersecting the last scattering surface could leave 
an imprint on the CMB ||, but these are unlikely to be 
observable with planned CMB probes Q] . 

If we drop the assumption of baryon symmetry, allow- 
ing for a lesser amount of antimatter than matter, then 
instead of a lower limit to the domain size, observations 
just place upper limits to the antimatter-matter ratio R 
at different distance scales. Indeed, it may be possible 
to have a small fraction R < 10~ 6 of antimatter stars in 
our galaxy ||. 

No antinuclei (with \ Z\ > 1) have ever been observed in 
cosmic rays. The Alpha Magnetic Spectrometer (AMS) 
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H to be placed on the International Space Station will 
look for antinuclei in cosmic rays, and if none are found, 
will place a tight upper limit on the antimatter fraction 
of cosmic ray sources. The AMS precursor flight on the 
Space Shuttle observed 2.86 x 10 6 helium nuclei but no 
antihclium j?J, giving an upper limit He/He < 1.1 x 10 _h 
on the antihelium-helium flux ratio in cosmic rays. 



B. Inhomogeneous baryogenesis 

Early work on antimatter regions in the universe (see 
the review by Steigman [Q) considered them as an ini- 
tial condition for the universe or tried to form them 
by separating matter from antimatter at a later stage 
Q . Later work is related to scenarios for inhomogeneous 
baryogenesis. 

There are many proposed mechanisms for baryogen- 
esis, including grand unified theory (GUT) baryogene- 
sis, electroweak baryogenesis, and Affleck-Dine baryo- 
genesis. The simplest versions produce a homogeneous 
baryoasymmetry, but simple modifications lead to an in- 
homogeneous baryogenesis which produces matter and 
antimatter regions 10 li) . See, e.g., the reviews by Dol- 
gov @. 

Inhomogeneous baryogenesis without inflation leads to 
a matter-antimatter domain structure with a very small 
distance scale. Models connected to inflation can lead to 
arbitrarily large distance scales. Some scenarios for GUT 
baryogenesis lead to an unacceptable large domain wall 
energy between the matter and antimatter domains, but 
other scenarios avoid this problem [ |ll| . 

Most studies of inhomogeneous baryogenesis have been 
for a globally baryon symmetric universe. As it has been 
recently shown || that the distance scale in this case 
would have to be at least comparable to the present hori- 
zon, the attention has shifted to models where the observ- 
able universe is baryon asymmetric, but could contain a 
smaller amount of antimatter [l^Jl7[| . 



C. Antimatter regions in the early universe 

On scales smaller than about 1 kpc [jl8|, antimatter 
regions would have annihilated by now but could have 
left an observable signature in the CDG spectrum, in the 
CMB spectrum, or in the yields of light elements from 
big bang nucleosynthesis (BBN). 

The smaller the size of the antimatter regions, the ear- 
lier they annihilate. Domains smaller than 100 matT = 
1 MeV, corresponding to a comoving (present) scale of 
6 x 10 8 km or 0.02 mpc, would annihilate well before 
nucleosynthesis and would leave no observable remnant. 

The energy released in antimatter annihilation ther- 
malizes with the ambient plasma and the background 
radiation, if the energy release occurs at T > 1 keV. 



If the annihilation occurs later, Compton scattering be- 
tween electrons heated by the annihilation and the back- 
ground photons transfers energy to the microwave back- 
ground, but is not able to thermalize this energy (because 
Compton scattering conserves photon number). The lack 
of observed distortion in the CMB spectrum constrains 
the energy release occurring after T = 1 keV to below 
6 x 10 -5 of the CMB energy |l9). This leads to progres- 
sively stronger constraints on the amount of antimatter 
annihilating at later times, as the ratio of matter and 
CMB energy density is getting smaller. Above T ~ 0.2 
eV the baryonic matter energy density is smaller than 
the CMB energy density, so the limits on the antimatter 
fraction annihilating then are weaker than 6 x 10~ 5 . 

For scales larger than about 100 pc (or 7 x 10 11 m at 1 
keV) the tightest constraints on the amount of antimatter 
come from the CMB spectral distortion, and from the 
CDG spectrum for even larger scales ||. 

We consider here intermediate distance scales, where 
most of the annihilation occurs shortly before or during 
nucleosynthesis, or after nucleosynthesis but before re- 
combination, at temperatures between 1 MeV and 1 eV. 
The strongest constraints on the amount of antimatter 
at these distance scales will come from big bang nucle- 
osynthesis affected by the annihilation process. 



D. BBN with antimatter 



Much of the early work on BBN with antimatter 



|T],^0 22 1 was either in the context of a baryon symmetric 
universe [ pb"f or for a homogeneous injection of antimatter 
through some decay process p2[ . 

Rehm and Jedamzik |^3| studied small antimatter re- 
gions, which annihilate before nucleosynthesis, at tem- 
peratures T > 80 keV. Because of faster diffusion of 
neutrons and antineutrons (as compared to protons and 
antiprotons) , annihilation reduces the net neutron num- 
ber, leading to underproduction of 4 He 0. This sets a 
limit R < 10~ 2 for the amount of antimatter in regions 
of size r A -lcmatT = 100 GeV (2 km at T = 1 MeV 
or 3 x 10 6 m at T = 1 keV). Our results for these small 
scales agree with ]2j| . 

We consider also larger antimatter regions, which an- 
nihilate mainly during or after nucleosynthesis. We 
have done detailed inhomogeneous nucleosynthesis cal- 
culations, where diffusion, annihilation, and nucleosyn- 
thesis all happen simultaneously. 

The case where annihilation occurs after nucleosyn- 
thesis was considered in f2l]| . Because annihilation of 
antiprotons on helium would produce D and 3 He it was 
estimated that the observed abundances of these isotopes 
place a comparable upper limit to the amount of anti- 
matter annihilated after nucleosynthesis. As we explain 
below, the situation is rather more complicated. 

We reported our first results in , where we had not 
included some features whose effect wc estimated to be 
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small. These included (1) antinucleosynthesis in the an- 
timatter region, (2) photodisintegration of other isotopes 
than 4 He, and (3) the dependence of the electromagnetic 
cascade spectrum on the initial photon spectrum from 
annihilation. We have now made the following changes 
to our computer code to take these effects into account. 

(1) We have added all antinuclei up to A = 4 and their 
antinucleosynthesis. Annihilation of these antinuclei pro- 
duce energetic antimatter fragments which may pene- 
trate deep into the matter region and annihilate there. 
Thus annihilation reactions occur also far away from the 
matter-antimatter boundary. 

(2) We have added photodisintegration of the lighter 
nuclei, D, 3 H, and 3 He. 

(3) We treat photodisintegration in more detail, espe- 
cially at lower temperatures, where use of the standard 
cascade spectrum is no longer appropriate. 

Below, all distance scales given in meters will refer to 
comoving distance at T = 1 keV. One meter at T = 
1 keV corresponds to 4.24 x 10 6 m or 1.37 x 10~ 10 pc 
today. Rehm and Jedamzik give their distance scales 
at T = 100 GeV. Our distances are thus larger by a factor 
3.0 x 10 s . We use h = c = k B = 1 units. 

The physics of the annihilation of antimatter regions in 
the early universe is discussed in Sec. II. We describe our 
numerical implementation in Sec. Ill and give the results 
in Sec. IV. We summarize our conclusions in Sec. V. 



II. ANNIHILATION OF ANTIMATTER 
DOMAINS 

A. Mixing of matter and antimatter 

Consider the evolution of an antimatter region, with 
radius ta, surrounded by a larger region of matter. We 
are interested in the period in the early universe when 
the temperature was between 1 MeV and 1 eV (age of 
the universe between 1 s and 30000 years). The uni- 
verse is radiation dominated during this period. At first 
matter and antimatter are in the form of nucleons and 
antinucleons, after nucleosynthesis in the form of ions 
and anti-ions. Matter and antimatter are mixed by dif- 
fusion at the boundary and annihilated. Thus there will 
be a narrow annihilation zone separating the matter and 
antimatter regions. 

Before nucleosynthesis the mixing of matter and anti- 
matter occurs mainly through neutron/antineutron dif- 
fusion, since neutrons diffuse much faster than protons. 
If the radius of the antimatter region is less than r ~ 10 7 
m, all antimatter annihilates before nucleosynthesis. In 
nucleosynthesis the remaining free neutrons go into He 
nuclei. The mixing of matter and antimatter practically 
stops until the density has decreased enough for ion dif- 
fusion to become effective at T s» 3 keV. 

Thus there are two stages of annihilation, the first one 
before nucleosynthesis, at T > 70 keV, the second well 



after nucleosynthesis, at T < 3 keV. The physics during 
the two regimes is quite different. The first regime was 
discussed in p3fl . We concentrate on the second regime 
in the following discussion. 

Hydrodynamic expansion becomes important at T ~ 
30 keV. At that time the annihilation of thermal electron- 
positron pairs becomes practically complete and the pho- 
ton mean free path increases rapidly. When the mean free 
path becomes larger than the distance scale of the baryon 
inhomogeneity, the baryons stop feeling the pressure of 
the photons, which had balanced the pressure of baryons 
and electrons. The pressure gradient then drives the fluid 
into motion towards the annihilation zone p5[p6| . This 
flow is resisted by Thomson drag. The fluid reaches a 
terminal velocity 26 
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Here e 7 is the energy density of photons, ot = 0.665 x 
10 _28 m 2 is the Thomson cross section, P is the pressure 
of baryons and electrons, and n* = n e - — n e + is the net 
electron density. With P « (ns + n e )T and |n*| « ub, 
we get a diffusive equation 
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for the baryon density, with an effective baryon diffusion 
constant due to hydrodynamic expansion 
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The hydrodynamic expansion alone does not cause 
mixing, but it significantly speeds up annihilation by 
bringing material towards the annihilation zone. The 
annihilation zone is surrounded by a depletion zone [Q, 
where the density of (anti)matter has decreased due to 
matter flow into the annihilation zone. The resulting 
pressure gradient maintains this flow. 

Antinucleosynthesis in the antimatter region produces 
antinuclei. The yields of these anti-isotopes are not in- 
teresting in themselves, since they are eventually annihi- 
lated. The annihilation of these antinuclei with nucleons 
produces energetic antinucleons and lighter antinuclei, 
which may penetrate deep into the matter region before 
annihilating. Thus, in addition to "primary" annihila- 
tion in the annihilation zone, there is also "secondary" 
annihilation outside this zone. 



B. Annihilation reactions 

The primary annihilation reactions occur at low ener- 
gies where reaction data is scarce or non-existent. The- 
oretically, the annihilation cross section is known to be- 
have as 1/v, when one or both of the annihilating parti- 
cles are neutral, and as 1/v 2 when both are charged. 
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More precisely, the theoretical nA cross section is 



^4^|^ Im( - a - ) -2Im 2 (- a , 



(5) 



where a s is the scattering length, q = /xv, \x is the reduced 
mass, and v is the relative velocity of the annihilating 
particles. 

The analogous expression for S-wave pA annihilation 

is MT 
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where ?y = —1/qB is the dimensionless Coulomb param- 
eter, B = l/Zfj,a is the Bohr radius of the antiparticle- 
particle system, a sc is the Coulomb-corrected scattering 
length and 
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There is laboratory data only for some of the rele- 
vant scattering lengths, and the uncertainties are large 
|50|,[n]]. From atomic data p^ , p3|p8[ , the pp system has 
Im(— a sc ) = 0.71 ± 0.05 fm. Recent experimental data by 
the OBELIX group |1 gives Im(-a sc ) = 0.62 ± 0.02 ± 
0.04 fm for pD and Im(-a sc ) = 0.36±0.03lo;ii for p 4 He. 

Primary annihilation is not sensitive to the annihila- 
tion cross sections, since annihilation is complete in the 
annihilation zone anyway. In secondary annihilation the 
A-dependence of the annihilation cross section is impor- 
tant, since it determines whether antinucleons annihilate 
with protons, which leads to no nuclear yields, or with 
4 He, producing D and 3 He. 

The yields of the annihilation reactions are important. 
Fortunately there is data on the most important reac- 
tion, antiprotons on helium |3q|, and also on some other 
reactions with antiprotons fl36p7| . 

The annihilation reaction between an antinucleon and 
a nucleus can be thought of as an annihilation of one of 
the nucleons in the nucleus. According to experimental 
data, an antiproton is twice as likely to annihilate on a 
proton than on a neutron in the nucleus |3^j39[] . 

The annihilation of a nucleon and an antinucleon pro- 
duces a number of pions, on average 5-6 with a third of 
them neutral . The charged pions decay into muons 
and neutrinos, the muons into electrons and neutrinos. 
The neutral pions decay into two photons. About half 
of the annihilation energy, 1880 MeV, is carried away by 
the neutrinos, one third by the photons, and one sixth 
by electrons and positrons. 

When an antinucleon annihilates on a nucleus, some 
of the produced pions may knock out some of the other 
nucleons, or in the case of larger nuclei, small fragments 



(p, D, 3 H, 3 He, 4 He). Some of the annihilation energy 
will go into the kinetic energy of these particles and the 
recoil energy of the residual nucleus. Typical energies are 
of order - 10 MeV. 

According to Balestra et al. |35| the average yields 
of low-energy p 4 He annihilation are 0.210±0.009 3 He, 
0.437±0.032 3 H, 0.07-0.19 D. This leaves about 0.7-0.9 
nucleons. 

Experimental data on the energy spectra of these emit- 
ted nucleons and fragments can be approximated by 
Ce~ E ' E °, where the average energy Eq decreases with 
the mass of the emitted particle. However the corre- 
sponding momentum is close to 350 MeV/c independent 
of mass Q. The momenta of the residual nuclei are 
smaller. Balestra et al. |35[] report a measurement on the 
momentum distribution of 3 He from p 4 He annihilation, 
with the mean energy corresponding to a momentum of 
198 MeV/c. Because of their large momenta these reac- 
tion products get spread over a large area, many of them 
escaping the annihilation zone, at least for a while. 



C. Thermalization of annihilation products 

The annihilation products lose their kinetic energy 
through collisions in the ambient plasma. Ions lose en- 
ergy by Coulomb scattering on electrons and ions and 
by Thomson scattering on photons. If the velocity of the 
ion is greater than thermal electron velocities, the energy 
loss is mainly due to electrons. At low energies scattering 
on ions becomes important. 

For an ion with E 3> T, we find that the energy loss 
per unit distance due to Coulomb collisions is 
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Here M, Z, and E are the mass, charge and energy of 
the incoming ion, T is the temperature of the plasma, m, 
z, and n are the mass, charge, and number density of the 
plasma particles, and A ~ 15 is the Coulomb logarithm. 
We assumed here that both the incoming ion and the 
plasma particles are non-relativistic. 

When the ion velocity is large compared to the thermal 
velocities of plasma particles, Eq. (||) simplifies into |h]] 
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and in the opposite case, [mE / {MT) <g; 1], into |^l[] 
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The energy loss in a plasma consisting of electrons and 
nuclei is thus 
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where we assumed m e /M <C 1, and used approximation 
(^) for scattering on nuclei. The index i labels different 
species of nuclei in the plasma. 

We plot the penetration distance d of 3 He and 3 H ions 
in a homogeneous plasma in Fig.|J as a function of E/T, 
in the absence of thermal electron-positron pairs. At 
large energies d cx E 2 , which corresponds to approxima- 
tion (^). We show also the effect of ignoring scattering 
on nuclei. 

The drag force exerted on the ion by photons is 



dE 4 / //., 
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The effect of this is negligible compared to Coulomb scat- 
tering. 
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FIG. 1. The penetration distance of a 3 He and a 3 H ion 
in matter with constant baryon density, i) = 6x 10~ 10 . The 
distance d is given in comoving units at T = 1 keV. The solid 
line is for a case where all baryons are in form of protons, 
and the dashed line for a 4 He mass fraction of 0.25. The 
dot-dashed line shows the effect of ignoring scattering on nu- 
clei. The penetration distance for an ion with charge Z and 
mass number A is obtained approximately by scaling the 3 H 
curve by a factor A/(3Z 2 ) vertically and by A/3 horizontally. 
For scattering on electrons (dot-dashed line) this scaling rule 
is exact. 

Neutrons lose energy through scattering on ions and 
electrons. Scattering on electrons is not important for 
T < 30 keV. The neutron loses a substantial part of its 
energy in each collision with an ion. The penetration 
distance is of order of the mean free path A = l/(an). 



Assuming -q — 6 x 10~ 10 , we find for neutron-proton scat- 
tering A ~ 4.7 x 10 9 m (77 keV)" 2 for a neutron with a 
typical 70 MeV energy. At T < 0.36 keV the mean free 
time of a 70 MeV neutron becomes larger than its life- 
time. The neutron is then likely to decay into a proton 
before thermalizing. 



D. Photodisintegration 

The high-energy photons and electrons from pion de- 
cay initiate electromagnetic cascades fl42|-[48|] . The dom- 
inant thermalization mechanisms for energetic photons 
and electrons are photon-photon pair production and in- 
verse Compton scattering 
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with the background photons jbb- The cascade proceeds 
rapidly until the photon energies E~ t are below the thresh- 
old for pair production, 



(14) 



where e 7 is the energy of the background photon. 

Because of the large number of background photons, a 
significant number of them have energies 3> T, and the 
photon-photon pair production is the dominant energy 
loss mechanism for cascade photons down to 45 
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Below this threshold energy, the dominant scattering pro- 
cess is photon-photon scattering down to H] 



E r = 
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Below this energy the dominating energy loss processes 
for photons are pair production on nuclei and Compton 
scattering on electrons. Inverse Compton scattering is 
still the dominant energy loss mechanism for electrons. 

When energy is released in the form of photons and 
electrons with energies well above E mlLX , the energy is 
rapidly converted into a cascade photon spectrum, which 
depends only on the total energy Eq injected, and is well 
approximated by [|45| , [47[ 
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The normalization factor is 
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Photon-photon pair production and scattering, and in- 
verse Compton scattering, are very rapid processes com- 
pared to interactions on matter, due to the large num- 
ber of photons. When the photon energies fall below E c 
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the mean interaction time rises drastically. The thermal- 
ization continues through Compton scattering and pair 
production in the field of a nucleus, in a time scale long 
compared with that of the cascade. The pair production 
cross section is [B9| 



3a ~ ,28, 2E 218, 

and the Compton cross section is (E 3> m e ) 
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Similarly, the energy of a photon from neutral pion de- 
cay (ir° — > 77) has a uniform distribution in the range 
^E n (l ± v n ). For a muon moving with velocity v and 
energy E^, the electron spectrum becomes 
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Photons with E > 19.9 MeV disintegrate 4 He, pro- 
ducing 3 He, and also D for E > 26.2 MeV. Above the 
energy E max the cascade proceeds so rapidly that photo- 
disintegration of nuclei is rare and can be ignored. The 
photodisintegration of 4 He begins at T — 0.6 keV, when 
E max becomes larger than the binding energy of 4 He. 
For T = 0.45-0.60 keV 4 He photodisintegration produces 
3 He (or 3 H) only, below T — 0.45 keV also D is produced, 
although with a smaller cross section. The photodisinte- 
gration of D begins earlier, at T = 5.3 keV, because of 
the smaller deuteron binding energy. The 3 He photodis- 
integration begins at T = 2.2 keV, 3 H at T = 1.9 keV, 
and 7 Li at T = 4.7 keV. 

During the second stage of annihilation, the mean free 
path of a photon at a given temperature is always larger 
than the distance scale of antimatter regions which an- 
nihilate at that temperature. We can therefore assume 
that the photons are uniformly distributed over space. 
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for0<£;< \E^{l-v). 

The electrons transfer their energy to background pho- 
tons through inverse Compton scattering. We calculate 
the scattering rate R for an electron with energy E c pass- 
ing through a thermal photon background, in the approx- 
imation E c 3> m c ^> T. Let E 1 be the energy transferred 
from the electron to a photon in one scattering. Using the 
Klein-Nishina cross section we get for a monochromatic 
photon background 
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E. Spectrum of annihilation photons and electrons where e = E 1 / {E c — £? 7 ), w = 4i? e 7 /m 2 , and e 7 is 

the energy of the background photons. Integration over 
the thermal photon spectrum gives the spectrum of up- 
scattered photons for one scattering, 



As the temperature falls the cascade spectrum moves 
to higher energies and, for T < 100 eV, it begins to over- 
lap the initial photon spectrum from annihilation. Then 
the lower part of this initial spectrum is no more con- 
verted to a cascade spectrum before photodisintegration, 
and the shape of the initial photon spectrum becomes 
important. 

In the pion's rest frame its direct decay products, pho- 
tons, muons, and muon neutrinos, have a single-valued 
energy, determined by conservation of energy and mo 
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mentum. The muon decays via fi 
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spectrum of the electron in the muon's rest frame is J50j] 

dn e E 2 ( 4E\ 1 . . , 

^ = 16^3 3-- , 0<E<-m, (cm.) (21) 



in the approximation m e /m A 



0. 



For the decay products of a moving pion, integration 
over directions yields an energy spectrum. The decay 
(tt^ — * /i* + u„) of a charged pion with velocity and 
total energy E v produces a muon with a uniform spec- 
trum in the range 



The average fractional energy loss (E 1 /E e ) in one scat- 
tering increases with increasing a. At a <C 1 the average 
energy transfer is (i? 7 ) = 3.60aE c . At large a the elec- 
tron loses most of its energy in one scattering. At the 
limit a< 1 the Klein-Nishina cross section reduces into 
the Thomson cross section. 

The electron scatters several times, losing a decreasing 
fraction of its energy in each collision. The process gen- 
erates a photon spectrum with most of the photons at 
low energies where dn/dEj oc Ej 3 ^ 2 . 
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In Fig. H we plot the spectra of electrons and photons 
from pion decay, for an exponential pion spectrum. We 
also show photon spectra resulting from inverse Compton 
scattering. 



100 




10 10 10 

E (MeV) 

FIG. 2. Initial spectra of photons and electrons from pion 
decay, for an exponential pion spectrum with the mean energy 
329 MeV. The dot-dashed line shows the photon spectrum 
from the decay of a neutral pion. The dashed line shows 
the electron spectrum from the decay of a charged pion. The 
electrons transfer their energy to background photons through 
inverse Compton scattering. The resulting photon spectra at 
temperatures 1 keV, 100 eV and 10 eV are shown by solid 
lines. 



F. Spallation of 4 He by energetic neutrons 

The average energy of a nucleon produced in p 4 He an- 
nihilation is ps 70 MeV. This is sufficient to disintegrate 
a He nucleus. Protons and ions slow down rapidly com- 
pared to the interaction time of nuclear reactions. Neu- 
trons thermalize much more slowly and may cause sig- 
nificant spallation of 4 Hc. 

Destruction of even a small fraction of 4 He may pro- 
duce 3 He or D in amounts comparable to the total abun- 
dance of these elements, but destruction of other ele- 
ments is significant only if a large fraction of the nuclei 
is destroyed. Thus only n 4 He spallation is important. 

For T < 100 eV the neutron mean time before spalla- 
tion becomes larger than the neutron lifetime and spal- 
lation gradually ceases. 



G. Lithium 

We do not expect any drastic effects on the 7 Li yield 
from antimatter regions. For small scales and large anti- 
matter fractions the reduction in the 4 He and 3 He yields 
cause an even steeper reduction in the 7 Li yield, but the 
4 He yield is a more sensitive constraint. 



For large scales, annihilation and photodisintegration 
of 7 Li is a small effect, just as for D and 3 He, compared 
to the large 3 He production from 4 He annihilation and 
photodisintegration. 

Since the standard BBN (SBBN) 6 Li yield is much 
below the 7 Li yield, 6 Li production from 7 Li annihilation, 
spallation, and photodisintegration could cause a large 
relative increase in the 6 Li yield. 

The 3 H and 3 He from photodisintegration and annihi- 
lation have large energies. They may react with 4 He to 
produce 6 Li and 7 Li before thermalizing. This nonther- 
mal nucleosynthesis may proceed via 3 H( 3 He) + 4 He — > 
6 Li+n(p), which has a threshold of 4.80 MeV (4.03 MeV) 
and is therefore not available for thermal nucleosynthe- 
sis, and it may result in a 6 Li yield much larger than in 
SBBN |5l|]. 



III. NUMERICAL IMPLEMENTATION 
A. General 

We use a spherically symmetric geometry where a 
spherical antimatter region is surrounded by a thick shell 
of matter. We assume equal initial densities rib — nj, i n 
both regions, such that the average net baryon density 
(jib) corresponds to (rj) = 6 x 10~ 10 . 

We give our results as a function of two parameters, the 
radius of the antimatter region ta, and the antimatter- 
matter ratio R. These parameters together with the net 
baryon density determine the initial local baryon density 
rib and the volume fraction fy covered by antimatter. 
The volume fraction depends only on the antimatter- 
matter ratio 



R = 



(1 - fy)n b 1 - /v 



(28) 



The initial baryon density rib is linked to the volume frac- 
tion through 

(n B ) = (1 - fy)n b - f v n b . (29) 

1/3 

The radius of our grid is L = ta / fy ■ We assume reflec- 
tive boundary conditions at the outer boundary of the 
matter shell. This models the situation where antimat- 
ter regions of radius ta are separated from each other by 
the distance 2L between their centers. 

For R <C 1, also fy <C 1 and ta <C L, so that we have a 
relatively small antimatter region surrounded by a much 
larger volume of matter. 

The annihilation creates a narrow depletion zone 
around the boundary between the matter and antimat- 
ter regions. An accurate treatment requires a dense grid 
spacing in this region. The position of the boundary 
moves with time. Therefore a fixed non-uniform grid is 
not adequate. We use a steeply non-uniform grid, which 
is updated at every time step. The number of grid cells 
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per unit distance is proportional to the gradient in baryon 
density. The total number of cells is kept constant. 

We include nucleosynthesis both in matter and anti- 
matter. In matter we follow the reactions up to A = 7, 
in antimatter up to A = 4. Our code includes 15 iso- 
topes: n, pJD, 3 H, 3 He, 4 He, 6 Li, 7 Li, 7 Be, n, p, D, 3 H, 
3 He, and 4 He. Heavier matter isotopes are included as 
sinks. 



B. Annihilation and diffusion 

Because of the large uncertainty or lack of data for 
most of the relevant annihilation reactions we simply use 



Annihilation reactions are represented by the differen- 
tial equation 



{crv) = cr 



(30) 



for the nn, np, np, and all nA, nA annihilation cross 
sections and 

(av) = C(v)a (31) 
for pp and all pA, pA, and AA annihilations. Here 



C(v) 



2Tr\ZiZ2\a/v 



1 -exp(-27r|ZiZ 2 |a/u)' 



(32) 



and we use oo = 40 mb. The velocity v in the Coulomb 
factor C{v) is the relative velocity, for which we use the 
thermal velocity v — y/3T/[i, where \x is the reduced 
mass of the annihilating pair. Wc also studied the effect 
of including an A 2 / 3 dependence in oq. 

We assume that nA and pA have the same nuclear 
yields, and that nA and pA have the corresponding an- 
tiyields. The most important pA reaction is p 4 He. For 
its yield we use 

p + 4 Hc -► 0.490n + 0.309p + 0.130D 

+ 0.437 3 H + 0.210 3 He, (33) 

3 H, 3 He yields from ||, 
|, and we assumed charge 



where we have taken the D. 
<r(pn)/o-(pp) = 0.42 from j| 
exchange has no net effect, to get the n and p yields. 
The nA, pA, nA, and pA yields for other nuclei than A = 
4 He are not important. We estimated yields for them by 
assuming that p(n) annihilation is twice as likely with p 
than with n in the nucleus p9"|]3§| ] , using the experimental 
p, D, and 3 H yields for p 6 Li and p 7 Li [ p7[ , and otherwise 
trying to mimic the p He data. 

There is no data on annihilation of an antinucleus on 
a nucleus. For simplicity we assume that the lighter nu- 
cleus is annihilated completely, and the remnants of the 
heavier nucleus go into He nuclei and nucleons, with 
equal number of protons and neutrons. Especially, anni- 
hilation of a nucleus on an antinucleus with equal mass 
number leads to total annihilation. 

Annihilation, nuclear reactions, and diffusion are 
solved together for better accuracy. Hydrodynamic ex- 
pansion, spreading of the annihilation yields, and photo- 
disintegration are treated as separate steps. We include 
diffusion of all ions and neutrons. 



dYk 
dt 



= - J2 (n B ) (oTv)Y k Yi = - £ G kl Y k Y u (34) 



where the indices k and I refer to the annihilating isotopes 
and Yfc = n k /{nB) is the relative abundance. 

We integrate Eq. ( |34"| ) over the time step At. We take 
the implicit equation 



(35) 



and linearize it into 

Y& ~ Y k = -^iSGuffiYl + Yfei - Y£Y?)AL (36) 



Here is the initial abundance and Y is the solution 
from the previous iteration step. 

We solve this equation iteratively by a modified 
Newton-Raphson method. The ordinary Newton- 
Raphson method (e.g., |5^| ) does not work in this case, 
since it often converges to an unphysical solution with 
negative Y. We stabilize the algorithm by introducing a 
parameter /i which initially is set to zero. We gradually 
increase the value of /i between iteration steps, until the 
solution has converged and /i = 1. 

In our code the hydrodynamic expansion is started at 
a constant temperature T = 30 keV. The results are in- 
sensitive to the starting temperature, because late times 
dominate the expansion. Equation (^) is solved for ns 
as an ordinary diffusion equation. The grid cells are then 
expanded so that the baryon distribution corresponds to 
the solution. 

Hydrodynamic expansion is not combined into the 
same matrix equation with diffusion and nuclear reac- 
tions, but is treated as a separate step. For this reason 
the convergence of the code requires a very small time 
step at late times, when the hydrodynamic expansion 
becomes important. This limited our ability to calculate 
with very large scales. For our largest scale ta = 10 11 m, 
we did not get a converged result for the CMB distortion, 
since it is sensitive to the lowest annihilation tempera- 
tures, although our results for the nuclear abundances 
did converge. 



C. Spreading of the annihilation products and their 
reactions 

We model the energy spectrum of a nucleus cre- 
ated in annihilation by an exponential distribution 
cxp(—E/E )/E - The spectrum is cut off at E = 10E . 
The mean kinetic energy Eo corresponds to momen- 
tum Po = 350 MeV/c for neutrons and protons, and to 
P a = 200 MeV/c for nuclei with A > 1. 
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Consider the spreading of nuclei produced during one 
time step, along a linear path. The spherical symmetry 
allows us to identify paths with same tangential distance 
ro from the symmetry center. Let F(E, s,ro)dro denote 
the cumulative spectrum of nuclei at distance s from the 
tangent point. The energy spectrum obeys the differen- 
tial equation 

dE 

F(E, s, ro)dro = F(E ^—ds, s — ds, r^jdr^ (37) 

ds 

+F (E)g(r)Airr 2 dr—. 

47T 

Here g(r) is the number of particles created per unit vol- 
ume at distance r from the center, Fq(E) is their initial 
spectrum, and solid angle fi(ro) = 27ryl — (ro/r) 2 picks 
directions which correspond to the tangential distance ro- 

We integrate Eq. ( p7| ) assuming a reflective boundary 
condition at the outer edge of the grid. Nuclei which fall 
below E m i n ~ T are considered thcrmalizcd. The formu- 
lae for the energy loss dE/ds due to various scattering 
processes were given in Sec. II. 

Ions lose energy through Coulomb scattering on elec- 
trons and ions, and Thomson scattering on photons. 
Neutrons lose energy through scattering on electrons and 
ions, or they decay and thermalize as protons. We include 
all these effects. Neutrons are allowed to scatter on an 
ion once, after which they are stopped. The strong and 
weak interaction neutron reactions included are listed in 
Table I. 



Reaction 



Ref. 



n+p total 
n+p total 
n+ 4 He total 
n+ 4 He -> S H + D 
n+ 4 Hc -f 3 H + p + n 
n+ 4 Hc -> D + p + 2n 
n+ 4 He -> 2D + n 
n+ 4 He -» 3 He + 2n 
n -> p 



1 55 
56 1 

57 (from inverse reaction), 

1 58 
1 58 
1 58 
58 

r = 886.7 s 



TABLE I. Neutron reactions and references to their cross 
section data. 



D. Nonthermal nuclear reactions 

We ignore spallation of nuclei by energetic nucleons for 
other nuclei than 4 He. Our results show that even 4 He 
spallation is a relatively small effect, which confirms that 
spallation of other nuclei can be safely ignored. 

We ignore in this work also the production of 6 Li by 
non-thermal 3 He( 3 H) + 4 He reactions 
incorporating it for future work [p9||. 
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but we are 



E. CMB distortion 

We calculate the ratio of injected energy to the CMB 
energy as 



W = 



2kcV 



f 



dp a 



Pcmb(T) dT 



'-dT 



(38) 



and require W < 6 x 10~ 5 to satisfy the CMB constraint. 
Here /?cmb is the energy density of the background radi- 
ation and yO ann is the energy density released in annihila- 
tion reactions in form of photons and electrons, averaged 
over space. Effectively, we are assuming complete ther- 
malization above T = 2 keV (redshift z ps 8.5 x 10 6 ) and 
no thermalization below it. We count into p ann half of the 
total annihilation energy. The other half disappears as 
neutrinos, and has no effect on nucleosynthesis or CMB. 



F. Photodisintegration 

We compute the initial spectra of electrons and pho- 
tons from pion decay following Sec. HE. We assume an 
exponential kinetic energy distribution for the pions, with 
mean total energy equal to 2m p /5.7 = 329 MeV. The 
electrons transfer their energy to background photons 
through inverse Compton scattering. We compute the 
spectrum of the upscattered photons using the Klein- 
Nishina cross section, assuming a thermal background 
spectrum and E e m e . We then redistribute the energy 
of the initial photons (upscattered and from ttq decay), 
whose energies are above E c into the standard cascade 
spectrum [Eq. (0)]. 

The photons in this resulting spectrum have then an 
opportunity to photodisintegrate. These photons may 
pair produce on a nucleus, Compton scatter, or photo- 
disintegrate nuclei. We allow an unlimited number of 
Compton scatterings for a single photon, but we remove 
the photon after the production of an e^ pair or a pho- 
todisintegration reaction. The created e^ pairs, as well 
as the background electrons which gain energy in Comp- 
ton scattering, will produce a second generation of non- 
thermal photons by inverse Compton scattering. These 
secondary photons are, however, much less energetic than 
the primary ones, and we ignore them. 

The photodisintegration reactions included in our code 
are listed in Table II. 

In pij we used the results of Protheroe, Stanev and 
Berezinsky (PSB) [[l6| for photodisintegration. PSB cal- 
culated the amount of 3 He and D produced per 1 GeV of 
energy released in the form of photons and electrons, as a 
function of redshift. However, their result does not apply 
for annihilation at low temperatures, when a significant 
part of the initial photon spectrum from annihilation is 
below the threshold for photon-photon pair production. 
In Fig. H we compare the PSB yields with the more de- 
tailed treatment described above which we are now using. 







0.15i 




T(eV) 

FIG. 3. Comparison of our photodisintegration 3 He and D 
yields (solid lines) with Protheroe et al. Pq| (dashed lines). 

We get less 3 He (D) than the PSB yield for T < 100 eV 
(50 eV). There are two reasons for this difference. 



Reaction Ref. 



D+7 -> 


p+n 


53 


(from inverse reaction) 


3 H + 7 - 


-> D+n 


60 


m 


3 H + 7 - 


-» p+2n 


60 




3 He + 7 


^D+p 


62 




3 He + 7 


-> 2p+n 


60 




4 He + 7 


-> 3 He+n 


63 





4 He + 7 -> 3 H+p same as 4 He + 7 -> 3 He+n 



4 He + 7 — > D+p+n |64| 

TABLE II. Photodisintegration reactions and references to 
their cross section data. 

First, our cross sections (Table II) differ somewhat 
from what PSB used. The main difference is for large 
photon energies E 1 > 200 MeV, where the 4 He photodis- 
integration cross section again becomes large and a pion 
is produced ("pion photoproduction" ) . PSB assumed 
large D and 3 He yields for these reactions. Available 
data |J5| gives very small cross sections for the 4 He+7 — > 
3 H + p, 4 Hc + 7 -> D + p + n, and 4 Hc + 7 -> D + D 
channels for E 1 > 200 MeV. Accordingly, we set the D 
and 3 He yields to zero in this range. Therefore we get 
lower D and 3 He production at low temperatures as the 
cascade moves to these higher energies. 

Second, for low temperatures the cascade energies 
move up to and beyond the energies of the initial an- 
nihilation photons. Since we only convert to the cascade 
those initial photons whose energy is above the cascade 
turnover E c , our photon spectrum for photodisintegra- 
tion does not move up further. Therefore the D and 
3 He yields become almost independent of temperature 
for T < 5 eV. In our nucleosynthesis runs most of the 
annihilation takes place for T > 5 eV, but for antimat- 
ter regions larger than ta ~ 10 11 m annihilation occurs 
at these lower temperatures and the photodisintegration 



contribution should become independent of r^. 

IV. RESULTS 

We show light element yields as a function of the radius 
of the antimatter region and the antimatter-matter 
ratio R in Figs. ||, || andg. We show also the temperature 
around which the annihilation is taking place. All results 
are for a net baryon density (77) = 6 x 10~ 10 . 

For scales smaller than ta = 10 5 m, annihilation hap- 
pens before the weak freeze-out, and has no effect on 
BBN. For scales between ta — 10 5 m and ta = 10 8 
m, neutron annihilation before 4 He formation leads to a 
reduction in the 4 He yield compared to standard BBN 
(SBBN). 

r A (P°) 
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FIG. 4. The yield of 4 He as a function of the antimat- 
ter-matter ratio R and the radius rA of the antimatter re- 
gions. The SBBN result, which is approached in the lower 
left corner, is Y p = 0.2484. The dotted lines show contours 
of T ann , the temperature at which half of the antimatter has 
annihilated. 

If annihilation is not complete before 4 He formation, 
it is delayed significantly because the neutrons have dis- 
appeared and ion diffusion is much slower than neutron 
diffusion. There will then be a second stage of annihila- 
tion well after nucleosynthesis, at T ~ 3 keV or below. 
This leads to a substantial increase in the yields of 3 He 
and D. 

Antimatter regions in the size range ta ~ 10 7 -10 8 m 
are annihilated in two stages. In the lower part of this 
range, practically all antineutrons diffuse out of the anti- 
matter region and are annihilated in the first stage, but 
neutrons diffusing towards the antimatter region manage 
to annihilate only an outer layer of the antiprotons before 
nucleosynthesis swallows the remaining neutrons. Thus 
the antimatter region that is left for the second stage of 
annihilation consists of antiprotons only. 
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FIG. 5. Same as Fig. @, but for 3 He. The SBBN yield is 
3 He/H = 1.06 x 10" 5 . 
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FIG. 6. Same as Fig. g, but for D. The SBBN yield is 
D/H = 2.70 x 10~ 5 . 



Larger antimatter regions, r^>4x 10 7 m, have also 
antineutrons left by the time of nucleosynthesis, and thus 
antinucleosynthesis, producing mainly 4 He, takes place 
in the antimatter region. The main significance of this 
is that p 4 He annihilation will later produce high-energy 
antinucleons, which penetrate deep into the matter re- 
gion before annihilating. Thus not all of the annihilation 
occurs in the annihilation zone ("primary" annihilation), 
but there is also a significant amount of "secondary" an- 
nihilation occurring in a large volume surrounding the 
annihilation zone. 

The main annihilation reaction during the second stage 
is p 4 He. It produces 3 He and a smaller amount of D. 
Because of their high energy, these annihilation products 
penetrate some distance away from the annihilation zone. 
Less than half of them end up in the antimatter region 
and are annihilated immediately. The rest end up in 
the matter region, but partly so close to the antimatter 
region that they are sucked into the annihilation zone and 
annihilated later (except for the largest scales studied). 

For ta ^ 5 x 10 7 m, part of the annihilation occurs be- 
low T = 0.6 keV where 4 He photodisintegration produces 
3 He and D. 

Thus there are two main contributions to 3 He and D 
production: annihilation and photodisintegration. We 
show these contributions separately in Figs. and ||. 

Figure 0a shows the net production of 3 He (including 
3 H) from all annihilation reactions. The most important 
3 He producing reaction is p 4 He. Another is n 4 He, where 
the antineutrons come from p 4 He annihilation. 3 He is 
destroyed primarily by p 3 He annihilation in the annihi- 
lation zone. 

Annihilation production of 3 He increases steeply from 
rA = 2x 10 7 m to 5 x 10 7 m as a larger part of the antimat- 
ter region survives till the second stage. For ta > 5 x 10 7 



m, the annihilation production of 3 He keeps increasing 
with scale, since the annihilation shifts to lower temper- 
atures where the 3 He produced in annihilation travels 
longer (comoving) distances, and is thus able to better 
survive annihilation. 

3 He produced in the matter region by the secondary 
annihilation is much more likely to survive and thus this 
secondary annihilation produces more, or at least a com- 
parable amount of, surviving 3 He than the primary anni- 
hilation in the range ta ~ 10 8 -5 x 10 9 m, where most of 
the 3 He from primary annihilation gets annihilated. In 
pij we did not include this secondary annihilation, and 
therefore we got a smaller annihilation contribution. 

Figure ^|b shows the net production of 3 He from all 
photodisintegration reactions. Photodisintegration of 
3 He sets in for > 2 x 10 7 m but has only a small 
effect. For > 5 x 10 7 m, photoproduction of 3 He from 
4 He overcomes 3 He photodestruction and increases up to 
r ~ 10 9 m, as a larger part of the cascade exceeds the 
4 He photodisintegration threshold. 

For ta ^ 10 10 m, photoproduction of 3 He decreases 
again, as the cascade keeps moving to higher energies. 
The photodisintegration cross sections for 3 He and D 
production are smaller at these higher energies, and be- 
cause the individual photons have higher energies there 
are fewer of them. For even larger scales, rA > 10 11 m, 
we would expect the photoproduction to stabilize as the 
cascade gets replaced by the initial annihilation spectrum 
(cf. Fig. § . 

The different dependence of these two contributions on 
T, and thus on rA, means that annihilation production 
dominates for r A = 2 x 10 7 -5 x 10 s m (T ann > 250 eV) 
and ta > 3 X 10 10 m (T ann < 10 eV), but photoproduc- 
tion dominates in the intermediate range r^ = 5x 10 8 - 
3 x 10 10 m. 
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FIG. 7. Contribution to the 3 He yield from a) annihilation 
and b) photodisintegration. To the left of the thick line in (b) 
the contribution is negative, since 3 He photodisintegration 
dominates over photoproduction. 



In Fig. ga, the feature at R > 0.1, r A = 10 9 -10 10 m is 
due to annihilation of the photoproduced 3 He. 

For D (see Fig. ||) we observe the same effects, with 
some differences. Annihilation produces about 5 times 
more 3 He than D, but D penetrates farther from the an- 
nihilation zone and thus survives better. Therefore the D 
yield from annihilation is less dependent on ta-, as most 
of the D survives already for smaller scales. The ratio of 
the net annihilation production of 3 He and D is there- 
fore less than 5, and approaches this number only for the 
largest scales, where finally most of the 3 He also survives. 

Photodisintegration of D begins already at T = 
5.3 keV, so it occurs always when the second annihila- 
tion stage is reached. Photoproduction of D from 4 He 
can only begin at T = 0.45 keV. Also the D yield from 
4 He photodisintegration is less than a tenth of the 3 He 
yield. Therefore D photoproduction overcomes photodis- 
integration only for scales ta ^ 3 x 10 8 m. 
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FIG. 8. Same as Fig. fj], but for the D yield. 

The third significant mechanism for D and 3 He pro- 
duction caused by annihilation is spallation of 4 He by 
the high-energy neutrons from p 4 He annihilation. For 
the scales ta = 10 7 -10 8 m its D and 3 He yields are about 
10% of that by annihilation reactions. For larger scales 
its relative importance falls off, as neutrons decay into 
protons, which are then thermalizcd, before encounter- 
ing a 4 He nucleus. 

Because of the large uncertainty about the annihilation 
cross sections in reactions involving other nuclei than just 
nucleons, we studied the effect of including an A 2 / 3 de- 
pendence in the cross section. This did not have a signifi- 
cant effect on the primary annihilation in the annihilation 
zone, but increasing the n He cross section increased the 
probability of secondary antineutrons annihilating 4 He 
instead of protons. Thus we got an increased 3 He yield 
for distance scales ta ~ 10 8 -5 x 10 9 m. Reducing the 
n 4 He cross section would have an opposite effect. 

Comparing our calculated yields to the observed abun- 
dances and the primordial abundances derived from them 
p6| , |67t , we obtain upper limits to the amount of antimat- 
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tcr in the early universe. We plot the limits from BBN 
and CMB on the antimatter-matter ratio R as a function 
of the radius of the antimatter region in Fig. |[ 

For small antimatter regions the limit comes from un- 
derproduction of 4 He. Using Y p = 0.22 as our lower 
limit to the primordial 4 He mass fraction, we obtain an 
upper limit R < 0.02-0.04 for r A = 0.6-20 x 10 6 m. 
Because this result is obtained from a calculation with 
the net baryon density (rj) = 6 x 10~ 10 , corresponding 
to the SBBN yield Y p = 0.248, a better way to state 
our 4 He constraint is that we allow a maximum reduc- 
tion of AY p = 0.028 from the SBBN result. Different 
assumptions on 77 and observed Y p could give a smaller 
acceptable AY p and thus a tighter limit on R. But this 
does not work in the other direction, since the 4 He yield 
falls very rapidly with increasing R. Thus the limit on 
R can hardly be relaxed from our stated value by using 
different observational constraints. 

At larger scales, ta > 2 x 10 7 m, the limit is set by 
overproduction of 3 He. There has been much uncertainty 
in the estimated primordial 3 He abundance, because of a 
large scatter in its observed abundances and uncertainties 
about its chemical evolution [ |66|j68| . Current knowledge 
suggests a probable primordial abundance of 3 He/H ~ 
10 -5 , with three times this value a reasonable upper limit 
Ipij . Thus we have used the constraint 3 He/H < 10~ 4,5 . 

The upper limit to R from 3 He falls rapidly as the 
distance scale is increased from 2 x 10 7 m to 10 9 m, where 
the limit becomes R < 2 x 10~ 4 . For even larger scales 
the limit is slightly relaxed but stays below 3 x 10~ 4 . 

Fig. ^ can be compared to Fig. 2 of Rehm and Jedamzik 
H or to Fig. 2 of Hi . Our 4 He yield is slightly larger 
and the corresponding limit to R weaker than in p3[ , 
because our net baryon density, (r/) = 6 x 10 -10 , is larger 
than the one used in (77) = 3.43 x 10~ 10 . Near 

rA ~ 10 8 m we now get a tighter limit on R due to a 
higher 3 He yield than we gave in (24). This is due to 
3 He production by secondary annihilation in the matter 
region, which was ignored in f24j| . 

These limits are stronger than those from the CMB 
spectrum distortion for scales r A < 10 11 m. We did not 
calculate the yields for larger scales, but the 3 He and D 
yields should become roughly independent of r^, since 
for these larger scales the primary annihilation products 
penetrate far enough from the annihilation region to sur- 
vive, and the spectrum responsible for photodisintegra- 
tion is the initial annihilation spectrum, so the depen- 
dence on the annihilation temperature disappears. The 
CMB limit should then become stronger than the 3 He 
constraint near the scale ta ~ 10 12 m. 



V. CONCLUSIONS 

We have studied the effect of antimatter regions of a 
comoving size ta ~ 10~ 5 -10 pc on big bang nucleosyn- 
thesis. Smaller antimatter regions annihilate before weak 
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FIG. 9. Upper limits from BBN and CMB to the antimat- 
ter-matter ratio R as a function of the radius ta of the an- 
timatter regions. The area above the solid lines is excluded 
by 4 He underproduction (Y p < 0.22) or 3 He overproduction 
( 3 He/H > 10" 4 - 5 ). The dashed line gives an alternative limit 
from using 3 He/D > 1 as the criterion for 3 He overproduction. 
The dot-dashed line is the limit from CMB distortion. 



freeze-out and are not likely to lead to observable con- 
sequences. Larger regions annihilate close to, or after 
recombination, and the amount of antimatter in such re- 
gions is tightly constrained by the CMB and CDG spec- 
tra. 

Regions smaller than r^~2x 10~ 3 pc annihilate be- 
fore nucleosynthesis. The annihilation occurs due to neu- 
tron and antineutron diffusion and leads to a reduction 
in the n/p ratio and thus to a reduction in Y p . Requiring 
Y p > 0.22, we obtain an upper limit R < few % for the 
primordial antimatter-matter ratio for antimatter regions 
in the size range ta ~ (0.1-2) x 10~ 3 pc. 

If the annihilation is not complete by nucleosynthesis, 
at T ~ 80 keV, it is significantly delayed, since all neu- 
trons and antineutrons are incorporated into 4 He and 
4 He, and (anti)proton diffusion is much slower. There 
will be a second stage of annihilation at T < 3 keV, 
when proton and ion diffusion finally become effective in 
mixing the remaining antimatter with matter. 

This second stage of annihilation leads to production 
of a large amount of 3 He and a smaller amount of D, 
through several mechanisms. 

Annihilation of 4 He with antiprotons in the annihila- 
tion zone separating the matter and antimatter region 
produces 3 He and D which are deposited some distance 
away from the annihilation zone. A large fraction of these 
annihilation products gets however sucked into the an- 
nihilation zone later and is thus annihilated. The sur- 
viving fraction increases with increasing distance scale, 
since this corresponds to a decreasing annihilation tem- 
perature. At lower temperatures the energetic ions from 
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annihilation penetrate larger (comoving) distances into 
the matter region. 

Antinucleosynthesis in the antimatter region pro- 
duces 4 He, whose annihilation produces antinucleons and 
smaller antinuclei. Of these, especially the antincutrons 
penetrate deep into the matter region, where they can 
annihilate 4 He producing 3 He, which has now a much 
better chance to survive. 

An important source of 3 He and D is photodisinte- 
gration of 4 He by the annihilation radiation. The large 
energy part of the initial radiation spectrum is converted 
into an electromagnetic cascade spectrum. The large- 
energy cut-off of the cascade exceeds the 4 He photodis- 
integration threshold when the temperature has fallen 
below 0.6 keV. For rvi > 3 x 10~ 2 pc most of the an- 
nihilation occurs below this temperature and thus the 
photoproduction of 3 He and D becomes important. Pho- 
todisintegration of 4 He is the dominant source of 3 He 
for rA ~ 0.1-10 pc. For larger distance scales the anni- 
hilation mainly occurs at lower temperatures where the 
photon spectrum shifts to higher energies where it causes 
less photodisintegration. 

Another source of 3 He and D is the spallation of 4 He by 
high-energy neutrons from annihilation reactions. This 
effect is at least an order of magnitude smaller than the 
ones discussed above. 

3 He and D are also destroyed by photodisintegration, 
but since the total amount of these isotopes is much less 
than that of 4 He, this is a small effect. 

For scales larger than ta ~ 2 x 10 -3 pc the tightest 
constraint on the primordial amount of antimatter is due 
to 3 He overproduction. p 4 He annihilation produces sev- 
eral times more 3 He than D and 4 He photodisintegration 
produces over 10 times more 3 He than D. For scales larger 
than ta ~ 0.1 pc, the requirement 3 He/H < 10~ 45 gives 
an upper limit R < 3 x 10~ 4 . 

Rehm and Jedamzik [[S9| have studied this same prob- 
lem and obtained results that seem to be in qualitative 
agreement with ours, but they find a lower 3 He yield. 
Their upper limit to R from 3 He overproduction is weaker 
than ours by about a factor of 2. They also criticize our 
use of 3 He/H as a constraint. Therefore we show in Fig. ^ 
also the constraint 3 He/D < 1, which is observationally 
more secure p8j. As can be seen from the figure, the 
limits to R stay essentially the same. By assuming that 
the low 6 Li/H observed in some Population II and disk 
stars is an upper limit to its primordial value, they ob- 
tain an even tighter limit on R from 6 Li overproduction 
J69| | . However, 6 Li is very fragile and is thus likely to be 
depleted in these stars. The main source of 6 Li is spal- 
lation by cosmic rays in the interstellar medium. Thus 
the primordial abundance of 6 Li based on observations is 
very uncertain, as noted also in (6{| , and could be much 
lower or much higher than the one observed. 

In conclusion, we have established nucleosynthesis con- 
straints on the amount of antimatter in the early universe 
which arc tighter, by a large factor, than those from the 
CMB spectrum, or any other known observational con- 



straint, for antimatter regions smaller than ~ 10 pc. 
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